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Abstract 



We comment on two recent calculations of the second order perturbative cor- 
rections in the heavy flavor semileptonic transitions within the Brodsky-Lepage- 
Mackenzie approach. It is pointed out that the results do not show significant 
enhancement either in the inclusive h — > c decays or in the exclusive amplitudes at 
zero recoil provided that the expansion parameter is chosen in a way appropriate to 
the kinematics at hand. The values of the second-order coefficients inferred from the 
BLM-type calculations appear to be of order unity in the both cases. Thus, in both 
cases no significant uncertainty in extracting V c b can be attributed to perturbative 
effects. The theoretical accuracy is mostly determined by the existing uncertainty 
in l/m 2 c nonperturbative corrections in the exclusive B — ► D* amplitude, and, to 
a lesser extent, by the uncertainty in the estimated value of the kinetic-energy ma- 
trix element /x^ in the case of r s i(6 — > c). The theoretical accuracy of the inclusive 
method of determination of V c h seemingly competes with and even exceeds the ex- 
perimental accuracy. 



* Permanent address 



In the last two years considerable progress has been achieved in the theory of 
nonperturbative effects in the heavy flavor decays. This theory has been applied to 
many problems from B physics, including such practically important application as 
precise determination of V c b from experimental data. Further work in this direction 
is under way. 

The current level of understanding of the nonperturbative effects is such that, 
surprisingly, calculation of the 'trivial' perturbative corrections due to hard gluon 
exchanges becomes the bottle neck of theoretical analysis. The perturbative correc- 
tions are often suspected to limit the theoretical accuracy. In most of the processes 
of interest only one- loop corrections are fully calculated at present. To match the 
level of accuracy obtained in the nonperturbative sector one needs to know two-loop 
0(a 2 ) terms. 

Complete calculation of these terms typically is a rather difficult task. Sooner or 
later it has to be addressed since without the complete two-loop results no analysis 
of the heavy flavor decays can be considered as fully conclusive. At present such 
calculations are not available, and one has to settle for less. 

It is known for a long time that b, the first coefficient in the Gell-Mann-Low func- 
tion, is a large numerical parameter in QCD. In many cases, when both 0(a s ) and 
0(a 2 ) corrections are known, the set of (a s /n) 2 terms is dominated by those graphs 
that are related to the running of a s . [] The sum of all these graphs is proportional 
to b(a s /n) 2 , while all other two-loop graphs yield (a s /n) 2 with a coefficient of order 
one. This observation gives rise to the so called BLM hypothesis || according to 
which a good idea of the (a s /7c) 2 terms can be obtained by calculating this "running 
a s " subset of the two-loop graphs while ignoring all other diagrams. 

In two recent stimulating papers the BLM approach has been used for estimating 
the a 2 corrections in the b — > c transition. In Ref. || (initiated by works |J) this 
correction was considered in the exclusive b —>■ c amplitude at zero recoil and was 
found to be small. (By small we mean that the coefficient in front of (a s /7c) 2 is 
of order unity.) A similar calculation was carried out in Ref. || for the inclusive 
width of the b — > c semileptonic decay with the conclusion that the a 2 terms are 
large. This result was interpreted (see ||) as a demonstration of an uncontrollable 
character of the perturbative series in the inclusive decays which allegedly blocks |5| 
accurate determination of V c b within the inclusive method. 

These two assertions above - small two-loop corrections in the exclusive case 
and large in the inclusive one - seemingly clash with each other if one combines 
them with the facts that (i) in the small velocity (SV) limit the total inclusive 
probability is completely saturated by the exclusive transitions B — > D* and B — > D; 
(ii) in the b —>■ c semileptonic decays we are actually rather close to the SV limit 
- in the essential part of the phase space the velocity of the D (D*) is small in 
the B rest frame. The proximity to the SV limit is substantiated by experimental 
data H according to which only 25 to 30% of the total probability leaks into non- 

1 Let us parenthetically note that in the Coulomb gauge these are merely the graphs with the 
polarization operator insertions in the gluon line EJ. 
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elastic channels (i.e. those other than B — > Dlv and B — > D*lv). The theoretical 
clarification of this experimental fact goes beyond the scope of the present paper 
and is discussed elsewhere. 

In this letter we argue that the paradox above is superficial. Analyzing the 
exclusive and inclusive channels in parallel and expressing the perturbative series 
in both cases in terms of one and the same coupling constant (the one which is 
physically relevant, see below) we eliminate the apparent qualitative contradiction 
between the results for inclusive and exclusive transitions and, moreover, demon- 
strate that the a 2 s terms are small, so that they cannot be reliably captured within 
the BLM approach per se. The present consideration, thus, gives further support to 
the observation || that the most accurate way of extracting V cb available now is the 
analysis of the inclusive semileptonic rate where the theory of the nonperturbative 
effects is much more predictive than in the exclusive decays. 

To elucidate our point we need to do some preparatory work. Let us start from 
the corrections in the exclusive b —>■ c amplitudes at zero recoil. The one loop result 
for the axial form factor tja at zero recoil is well known @, 

VA = l + — [ In ; 1 

7r \m b — m c m c 6 J 

The analysis of Ref. [ID] establishes that the normalization point of the gauge 
coupling a s that is physically adequate in this case is 

m = y/m c m b . (2) 

In other words, if one expresses the perturbative series for t]v,a i n terms of a s (mo) 
no large numbers appear in the expansion coefficients at order a 2 s . (Below, if not 
indicated explicitly, a s means, by definition, a s (mo).) This leads to a numerical 
estimate Q t\a — 0.965. 

At the two-loop level it was found in Ref. [[J that in the MS scheme 

VA = 1-0.431^- 1.211 (*)\ (3) 

The values of the coefficients above refer to the ratio of the quark masses 

z = m c /m b = 0.3 . (4) 

We will adopt this value of z throughout the paper for the sake of definiteness. All 
quantities referring to the MS scheme are marked by bars. The following generic 
notation will be used: 

a s I ' ol s \ 2 ,_ N 
r) A = l + a 1 ha 2 — + ... (5) 

7T \ 7T / 



2 Some partial higher order calculations of the b — > c zero recoil amplitudes which existed in the 
literature previously jllj |l2|, |(| and gave rise to the estimate tja = 0.986 ± 0.006 were shown [10 
to be irrelevant. 
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It is quite evident that the perturbative coefficients starting from the second 
order depend on the scheme used to define the strong coupling. The number for the 
two-loop coefficient in Eq. (|3|) is given in the MS scheme often used in calculations 
performed in dimensional regularization. This scheme uses certain rather ad hoc 
subtraction constants which are introduced for technical purposes only, and is quite 
unphysical. Moreover, such a choice is definitely unnatural for the calculations 
where one accounts only for the running of a s . The so-called V scheme of BLM 
is physically preferable here reflecting the real size of the effect. As a matter 
of fact, the actual calculations within the BLM approach are routinely performed 
just in this scheme (for details see |j~J|, [3|, Q), and only at the final stage is the 



series reexpressed in terms of a s in the MS scheme and d 2 introduced. It is more 
adequate to use the V^-scheme strong coupling, at least at the current stage when 
the calculations beyond the BLM level have not been done. In what follows we will 
work with the V^-scheme a s , unless otherwise indicated; a s no bar refers to the V 
scheme while a s is reserved for the MS coupling. 

The one loop relation between a s and a s is given by 




7T 

and, therefore, the perturbative expansion in the MS scheme 

t} A = 1 + ai h a 2 + ... 

7T V 7T / 



(6) 



takes the form 



7T 




VA = 1 • 77,^ ■ (a,- -^,); • -| + ... (7) 



so that ai = ai and a 2 = a 2 — j^bai . The first three terms in the expansion of t/a, 
Eq. (^|), then look as follows for m c /m h = 0.3 : 

VA = 1-0.431^5) + 0.286 f^V+ - (8) 



7T \ 7T 

where m is defined in Eq. (0). We see that the calculated second-order coefficient 
is in fact noticeably less than 1, i.e. clearly not enhanced. For this reason it cannot 
be taken even as an estimate of the second order result because other terms beyond 
the BLM approximation are expected to be of order 1. 

Now let us proceed to the inclusive semileptonic widths to bring the correspond- 
ing analysis in line with the exclusive case. 

Assume for a moment that (m& — m c )/(mb + m c ) 1; then the SV limit is 
obviously parametrically guaranteed, and the inclusive width is directly expressed 
in terms of tja and r)v, the zero recoil form factors of the vector and axial currents: 




(9) 
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This formula is nothing else than the text-book expression for the neutron /3-decay 
width (with the electron mass neglected). Moreover, in the SV limit the deviation 
of the zero recoil vector form factor rjv from unity is also proportional to the same 
smallness parameter having the meaning of the square of the typical velocity in 
the decay, and thus can be neglected as well. Therefore in the limit when (m& — 
m c ) I (mi, + m c ) <C 1 the perturbative expansions for the inclusive width and for the 
B — ► D* Iv zero recoil amplitude squared are directly related (essentially the same). 

In the real world m c is significantly smaller than m^. However, as was mentioned 
above, in the semileptonic b — > c decays one effectively is not far from the SV limit. 
The experimental evidence in favor of this fact is the approximate saturation of 
the experimental inclusive probability by the two elastic channels. It was pointed 
out in |T{J that the coefficients in the perturbative expansion of 77^ do not deviate 
drastically from their values at m c = mj; this observation was confirmed by explicit 
calculation in the BLM approximation in Ref. 0]. Since near the SV limit the 
elastic channel saturation is complete (this, of course, refers also to the perturbative 
corrections), it is advantageous to write the perturbative series for the inclusive 
width in terms of the same coupling constant as in the exclusive case, a s (mo). This 
choice rules out parametrically large two-loop coefficients in 77's; by the same token 
it will kill such coefficients in the inclusive width. 

Let us examine the result for the inclusive width || more carefully taking advan- 
tage of the proximity to the SV limit. Let us write down the perturbative expansion 
for the inclusive width in the form 




where Zq(z) is the well known phase space factor. Then one has 




Comparing the expansions for tja and for we get 

h(z) = ^(1) + O ((1 - zf) = -\ + 0({l-zf) , 

= ^(1) + + O ((1 - zf) = 3 -a 2 (l) + 1 + O ((1 - zf) . (12) 

It is essential that terms 0((1 — z)) are absent. Examination of Eq. (|T2|) shows 
that there are no reasons to expect too big a value of hi even for the actual masses 
of b and c. It is worth noting here that Eq. ( pT|) assumes the use of the "pole" 
masses for the b and c quarks, and the corresponding definition of the perturbative 
factor /tp in Eq. fllCf); otherwise Eq. (|9]) does not hold. The notion of the pole mass 
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cannot be consistently defined to all orders in the coupling constant, see Refs. |J. 
Perturbatively it can be defined to a given (finite) order, and for our limited technical 
purposes of discussing the two-loop perturbative corrections it is legitimate to use 
the corresponding "two-loop pole" masses. 

Of course, for Eq. flT2"D to hold the strong coupling must be defined in the same 
way in both cases, inclusive and exclusive. Taking into account the experimental 
proximity to the SV limit in the semileptonic b c decays, we will express the 
perturbative expansion for T s i in terms of a s (mo), which will allow one a direct 
comparison with the exclusive decays. 

Recently, calculation of the a 2 s effects in T s \(b c) associated with running of 
ot s has been done || based on a trick suggested in |13[]. For m c /m b = 0.3 it was 
obtained that 

4 = 1-1.67^- 15.1 f^V . (13) 

According to Eq. (||) and the one- loop evolution of a s (n) this expansion in terms of 
a s (mo) (the proper choice for the kinematics at hand) takes the form 

4 = l_l. 67 ^!^l_5.1l('M 2 + ... ~ 1 - 0.16 - 0.047 + ... (14) 



7T \ 7T 



where we used for numerical illustration the value a s (mo) = 0.3 corresponding to the 
choice a s (m ) = 0.24 adopted in Ref. Perturbative corrections, thus, do not show 
enhancement and, rather, seem to be well under control; no trace of uncontrollable 
10% second-order corrections claimed in is seen. We should note once again that 
the second-order coefficient in Eq. ([14]) is not very large; it becomes even smaller 
if one considers Kr instead of 4 which is reasonable for the purpose of comparison 
with 7] a (see below). If so, this result cannot be taken too literally because other 
second-order corrections that are beyond BLM are expected to be of a similar size. 

Now we turn to the question of theoretical uncertainties in determining 
Following the existing tradition to phrase the perturbative corrections to the B — > 
D* tv rate in terms of corrections to the amplitude rather than to the rate itself, 
we also write the perturbative expansion for the square root of the width, kt = 
^■pperiypo^i/2. form is directly related to extraction of \Va\: 

r] A = l- 0.431 a '( mo ) + Q.286 ■■■ = 1 - 0.041 + 0.0026+ ... 

7T V 7T / 

Kr = 1 - 0.835 a ^ m °^ -2.90 (— V+ ... = 1 - 0.08 - 0.026 + ... (15) 



7T \ 7T 



It is obvious that no significant theoretical uncertainty can be attributed to the 
perturbative corrections in both cases. Moreover, the genuine relation between 
the second-order coefficients can be established only after the complete two-loop 
calculations are made: the BLM-type calculation does not give a large contribution 
that would safely dominate the second-order result. 
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Continuing the analysis of the numerical aspect of the uncertainties it is ap- 
propriate to note that the nonperturbative corrections to the inclusive semileptonic 
widths are calculated |T3|, [T^, [RJ and give the correction 

(<M nonpert * -2.5% . (16) 

This is smaller than the first order perturbative correction but dominates the un- 
certainty in the a 2 s terms. A different situation occurs in the zero recoil B — > D* 
amplitude. The leading nonperturbative corrections 0(l/m^) are not fully known 
here (in a model-independent way); a reasonable estimate of these corrections is 
about —10% |S| with the uncertainty of about ±4%. Thus, they are definitely larger 
than even the first order radiative corrections and, therefore, the exact value of the 
second-order perturbative coefficient is not of much practical relevance here since it 
is not abnormally large. 

To complete our discussion of the purely perturbative uncertainties in b — > c 
transitions, let us mention that at present the value of the strong coupling at the 
scale 2.5 GeV is known, apparently, with the accuracy ~ 10%. The correspond- 
ing uncertainty in i]a is thus ~ 0.004 and is definitely irrelevant in view of large 
nonperturbative uncertainties. The corresponding uncertainty in Kr which enters 
determination of \V c b\ from the inclusive width constitutes about 0.007; being rather 
small and negligible at present, it may contribute a noticeable (through evidently 
not dominant) part in the overall theoretical uncertainty in the near future. At 
the moment the lack of the definite knowledge of the exact second-order coefficient 
also does not seem to limit the theoretical accuracy of this most precise method of 
determination of |K&|- 

Two brief remarks concerning the nonperturbative corrections which were men- 
tioned in passing are in order here. The theoretical basis for their calculation is 
provided by the Wilson operator product expansion (OPE) |17|]. In calculating 
r s i(6 — > c) OPE is used in the Minkowski domain which assumes duality. Deviations 
from duality (they die off exponentially) are controlled by the behavior (divergence) 
of the high-order terms in OPE [|TB||; they can not be found from purely theoret- 



ical considerations at present. On the basis of the phenomenological information 
available it is reasonable to believe that these deviations are negligible in the b —>■ c 
transition (unlike in c — > s decays ||19|| ). They are not included in the uncertainties 
in the nonperturbative terms which were quoted above. Another aspect of OPE 
ignored so far in this letter is the need for introducing a boundary point \x in any 
consistent OPE calculation with the subsequent separation of the short- and large- 
distance contributions. The latter must be subtracted from the coefficient functions 
calculated perturbatively since they are accounted for in the matrix elements of 
higher- dimension operators (condensate terms). In the calculations of the pertur- 



bative corrections to r/ A || and, especially, to the inclusive c — > s iv width [f20| the 
sizable part of the second-order coefficient comes from the low-momentum region 
of integration where perturbative expressions for gluon and quark propagators are 
not applicable. Therefore the numerical estimates derived from those calculations 
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have no direct physical meaning: the impact of this region is completely accounted 
for by nonperturbative contributions and must be excluded from the perturbative 
coefficients in the proper treatment (see |TJ, |2T|). If it is done, the size of perturbative 
corrections in c — > s Iv decays reduces essentially. In the case of t]a this modification 
is numerically larger than the whole second-order correction estimated in Ref. 0]: 



3n ' mi ' 3m c rrib 



*u(/i)^ + ^V \— 2 +— 2 +^^\ (17) 



where /i is the infrared cutoff (for more detail see [^T]]). The question of the consistent 



simultaneous account for both nonperturbative and perturbative effects [^TJ becomes 
important already at the level of the second-order perturbative corrections when the 
former are governed by the mass of c quark, due to the existing numerical hierarchy 
among the corrections in this case. Although the need for such consistent treatment 
based on the Wilson OPE is acknowledged, practically this aspect is usually ignored 



in HQET calculations (see, e.g., [TJ|, |22fl ) 



To summarize, we demonstrated that using the proper perturbative expansion 
parameter makes the second-order perturbative corrections equally small both in 
the b — > civ amplitude squared at zero recoil and in the inclusive b —> civ width. 
The existing BLM-type calculations do not show any particular enhancement of the 
second-order effects associated with the running of a s if the scale is chosen in a 
way appropriate for the SV kinematics. For this reason the recent results |3|, [| on 
the second-order perturbative corrections in the inclusive and exclusive heavy quark 
transitions can be viewed mainly as the indication that these corrections are not large 
and the perturbative series are well controlled in both cases. More precise results 
can be obtained only by complete two-loop computations; they are necessary for the 
ultimate improvement of the theoretical precision in the most accurate method based 
on the inclusive semileptonic width. For exclusive transitions this can not improve 
the accuracy of extracting V c b because the uncertainty is by far dominated by rather 
poorly known nonperturbative corrections. Still, the full two-loop calculation seems 
to be desirable in this well: one can use the value of V c b obtained from the 

inclusive widths to measure the exclusive form factor with sufficiently high accuracy 
and, thus, to obtain valuable information on the structure of nonperturbative effects 
in the b —>■ c transitions in the situation when they are significant. 

We emphasize that the concrete method of determining V c b from the inclusive 
width at present is based on the fact that the expression for the width valid in the SV 
kinematics works well for the actual quark mass ratio - the fact known theoretically 
for a long time and now substantiated by experimental studies of the final charm 
states showing the relevance of the SV limit. Therefore, the consistent analysis of 
corrections requires using the expansion parameters inherent to the SV kinematics. 

The result of recent calculations of potentially dominant second-order corrections 
H shows that the perturbative uncertainty affecting the extraction of V c b now does 
not exceed a percent level and is similar to the perturbative uncertainty in the 
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exclusive zero recoil transition, contrary to rather ad hoc claims existing in the 
literature |6|, |3|. This uncertainty is smaller than the impact of the nonperturbative 
corrections which are reliably calculated for the inclusive b — > c width. On the 
contrary, the nonperturbative corrections seem to be much larger in the exclusive 
decays and, moreover, they carry a significant uncertainty that can not be eliminated 
in a model- independent way at present. This limits the accuracy of the exclusive 
method of extracting V c b but, also, offers opportunities of learning more about details 
of the nonperturbative dynamics. 

It seems justified to say that with the recent progress in the theoretical treatment 
of the heavy flavor transitions the theoretical accuracy in determination of from 
T sl (B) achieved by year 1995 exceeds the existing experimental accuracy. 
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